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Mark 4, CORl/STEREO and LASCO/SOHO data analysis shows that 
at least a portion of type II radio bursts observed in the corona occurs in 
the presence of a CME, but in the absence of a shock ahead of them. A drift 
current instability in the CME frontal structure is discussed as a possible cause 
of such bursts. 
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1. Introduction 

Experiments have revealed the existence of a turbulence zone ( "foreshock" ) ahead of 
both the near-Earth shock [Cairns and Robinson, 1999] and interplanetary shocks trav- 
elling in the heliosphere [Bale et al., 1999]. In the foreshock, there are fluxes of energetic 
particles (electrons and ions) which move away from the front along undisturbed magnetic 
field. They are the most energetic portion of the plasma heated in the shock front. As 
a result of beam instability evolution, electron fluxes excite Langmuir oscillations at the 
plasma electronic frequency. Due to the Rayleigh and Raman scattering [Zaitsev, 1965], 
these oscillations are transformed, respectively, into the first and second harmonics of the 
type II radio emission at the single and double plasma electron frequencies [Kuncic et ai, 
2002]. Direct observations of the synchronous onset of increased Langmuir oscillations 
ahead of the interplanetary shock front as well as type II radio bursts at the same fre- 
quencies [Bale et al., 1999] confirm the existence of this process. It is suggested that a 
similar process may take place for shocks excited by coronal mass ejections (CMEs). No 
direct evidence has so far been obtained, however, due to difficulties involved in recording 
a shock front in the corona. Recently developed methods for directly recording shock 
fronts in the corona [Eselevich and Eselevich, 2008; Eselevich, 2010] allow one to come 
closer towards solving this problem. 

This study aims to demonstrate that, at least, a portion of type II radio bursts observed 
in the corona occur in the presence of CMEs, but in the absence of a shock ahead of them. 
A possible mechanism is discussed for the origin of such an emission associated with a 
drift current instabihty evolving in the leading edge of the CME frontal structure. 
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2. Input data and their representation 

Our analysis employed coronal images obtained by Mark 4 (Mauna Loa Solar Observa- 
tory littp://mlso.liao.ucar.edu), CORl (STEREO) [Howard et al, 2008] and LASCO C2 
(SOHO) coronagraphs [Brueckner et al, 1995]. The data were represented in the form of 
difference brightness (polarization or full) AP(i) — P{t) —P{to), where P{t) is the coronal 
brightness at moment t corresponding to the event under consideration, and P{to) is the 
undisturbed brightness at fixed moment to, selected as being well before the event start. 

Difference brightness images were used to inspect the dynamics of a CME and its 
disturbed zone. Representations in the form of isohnes AP, as well as sections both along 
the Sun radius at fixed position angles PA and non-radial sections at various instants 
were employed. Position angle PA is measured counterclockwise from the North Pole in 
the images. 

3. Data analysis method 

The analysis method is based on the results of recording the shock front in the corona 
[Eselevich, 2010]. Let us briefly review these results on the example of two CMEs with 
greatly different velocities: 5 May 1997 and 20 September 1997. Both these events are 
limb CMEs. Velocity with respect to the solar wind, u, for these CMEs is, respectively, 
150 km/s and 700 km/s {u — V — Vsw, where V is the measured CME velocity, and Vsw 
is the solar wind velocity in the streamer belt). Figures lA and B depict these CMEs in 
the form of difference brightness isohnes, at a certain instant for each event. The frontal 
structure is nearly circular in shape (dashed circumferences in Figures lA and IB) for 
both the CMEs. The center of the circle, O, is Rc away from the center of the Sun along 
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the position angle PAc- The direction from the frontal structure center is defined by 
angle a, measured counterclockwise from the CME direction (Figure IB). 

The differences between the slow and fast CMEs are conspicuous in both the images 
(Figures lA and B) and the difference brightness profiles AP{R) plotted along the T+ 
propagation direction (section 1 in Figures lA and B). These sections in Figures IC and 
D are plotted in the frontal structure coordinate system, which allows the formation of 
a shock discontinuity to be traced practically unambiguously in both space and time [M. 
Eselevich, 2010]. One can see from Figure 1 that: 

In the case of a slow CME (u ^ 150 km/sec), the difference brightness isolines 
are extend in the CME propagation direction. The difference brightness profile AP{K) 
in the disturbed zone ahead of the CME (black circles in Figure IC) smoothly decreases 
with distance. There is practically no disturbed zone in the transverse direction a = —70° 
(section 2, light circles). 

In the case of a fast CME (u « 700 km/sec), the difference brightness isohnes are 
nearly circular in shape. A discontinuity (jump) on the scale of ~ 0.25 R0 (where Rq is 
the solar radius) is observed in the frontal part of the disturbed zone in the AP{R) profile 
(black circles in Figure ID). Its position is indicated by a dashed bold hne in Figure IB. 

This discontinuity was identified with the shock front based on the analysis of more 
than thirty limb CMEs with velocities V ~ 200 -7- 3000 km/s. The analysis allowed the 
following steady regularity (law) to be discovered: "A disturbed zone extended along the 
CME propagation direction exists ahead of a coronal mass ejection when its velocity u 
with respect to the surrounding coronal plasma is below a certain critical velocity, uq- 
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Shock formation ahead of the CME frontal structure in the vicinity of CME propagation 
direction depends on whether local inequality u{R) > uc ^ V^(-R) holds, which can be 
true at various distances R> 1.5 from the Sun center. Here, Va{R) is the local Alfven 
velocity of slow solar wind in the streamer belt calculated in Mann et al. [1999]". 

The local Alfven velocity, Va{R), is shown as a sohd curve in Figure 2, with Vsw{R) 
plotted as a dotted curve. These statements serve as a basis for identifying shocks in this 
study and their relation to the presence or absence of type II radio emission. 

4. Event selection 

The following three CMEs were selected for the analysis: 25 March 2008, 31 December 
2007 and 3 June 2007. These events were investigated in detail in Gopalswamy et al. 
[2009] and have the following pecuharities: 

1. These CMEs occur close to the limb when observed by all the instruments ("limb" 
events). The velocity values measured by Mark4, CORla, b and C2 differ only shghtly 
and are close to the radial velocity. 

2. The white-hght CMEs have the simplest and three-part structure most convenient 
for investigation: frontal structure (FS), decreased brightness region (cavity) and bright 
core [Illing and Hundhausen, 1986]. 

3. Their propagation in the corona was accompanied by meter and decameter type II 
radio bursts. 
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5. Analysis of shock presence/absence for three CMEs associated with type 
II radio emission 

CME 1, 25 Mcirch 2008. Figure 3 shows the difference polarization brightness profiles 
AP{R) along the direction close to the CME motion axis {PA ^ 100°), drawn using 
the Mark 4 data for three consecutive instants. Light circle profiles correspond to 18:41, 
directly prior to the CME appearing within the Mark 4 field of view (undisturbed corona) . 

The position of the point from which it was drawn (parameters PAc and Rc specifying 
the approximate center of the frontal structure) and the direction (angle a) are indicated 
for each profile. For convenient comparison, the distances R from the Sun center are 
marked on the x axis. Additionally, the dashed bold line in the upper panel of Figure 3 
depicts the profile AP{R) at 18:47 in the direction a = 95° (sideward with respect to the 
CME motion direction). Its x position is calculated taking into account that the shape 
of the frontal structure is close to a circumference. In this direction, there is practically 
no disturbed zone ahead of the frontal structure, the frontal structure boundary (current 
sheet) having the size Si ~ 0.03 Rq (light gray on the profile), i.e. being close to the 
Mark 4 spatial resolution K 0.02 Rq. The brightness jump SPj at the frontal structure 
boundary has an amplitude close to the brightness jump amplitude SPf at the shock 
front. 

In the region of small a, there is a disturbed zone ahead of the frontal structure (shown 
by crosshatching in the profiles), bounded by a shock in the forefront (dark gray). (Shock 
propagation is seen in the following two plots in Figure 3). The disturbed zone being, in 
this direction, behind the shock front, it appears to include shock-heated plasma. 



DRAFT 



January 21, 2013, 3:29pin 



DRAFT 



X ^SELEVICH AND ESELEVICH: A FORMATION MECHANISM FOR THE TYPE II RADIO EMISSION 

The shock front velocity, u, exceeds Va (black squares in Figure 2) at all the distances. 
According to Gopalswamy et al. [2009] the meter type II radio emission is recorded when 
the forefront of the CME, corresponding to the shock, is at 1.7 R© < R < 2.2 Rq (this 
part is crosshatched and labeled 1 in Figure 2). The decameter type II emission ceases 
to be recorded when the shock reaches R ~ 3.6 (vertical straight line in Figure 2). 
Hence the type II radio burst is accompanied by CME driven piston shock propagation, 
in this case. 

CME 2, 31 December 2007. Figure 4A shows the difference brightness profiles based 
on the CORla data at consecutive instants for this event. Unhke the profiles in Figure 3, 
they are drawn in the frontal structure coordinate system, namely [Eselevich M, 2010]: 

1. each profile is normalized to its maximum value in the vicinity of the frontal structure 
(i.e. the maximum value of the normalized profile is 1); 

2. each profile (except the very earliest one) is shifted along the x axis so that the 
positions of the frontal structure coincide for all the profiles. 

Horizontal hatching shows the position of the frontal structure (the region of the nor- 
malized profile values > 0.5) in Figure 4A. For the initial instant 01:00, the difference 
brightness profiles for two angles a: 30° (diamonds) and 50° (dashed bold curve) - are 
shown in Figure 4A. The profiles are practically similar. No disturbed zone is perceptible 
ahead of the frontal structure for both of these directions. The steepest segment of the 
current sheet at the frontal structure boundary (light gray) in the direction a — 50° has a 
characteristic spatial scale Sj ~ 0.02-0.03 R0, comparable to the CORl spatial resolution 
K 0.016 R0. At the next instant 01:05 (asterisks in Figure 4L) the disturbed zone 
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increases only slightly in the direction a = 30°. Thus the current sheet of the frontal 
structure remains practically undisturbed prior to this instant. Its velocity u only slightly 
exceeds Va (light triangles in Figure 2). At the subsequent instants (01:10 and 01:15, not 
shown in Figure 4L), the disturbed zone continues to gradually increase, a jump - a shock 
with the front width Sp - eventually forming in the very forefront of the zone. It is shown 
in dark gray color for the instants 01:20 (crosses) and 01:30 (triangles). Such a change in 
the profile is absent in the direction a — 50°. 

Unlike the faster CME 1, where a horizontal segment of the AP{R) profile (i.e. 
AP{R) Ki const. Figure 3) is observed directly behind the front, this event exhibits no 
steady state directly behind the shock front even for the very latest instant 01:30 (trian- 
gles), since the difference brightness increases as the distance decreases. This may possibly 
mean that the shock front has been in the forming stage up to the instant in question. 
The encircled solid triangles correspond to it in the plot u{R) in Figure 2. It is only at 
large distances R > 4 Rq (non-encircled black triangles in Figure 2) that a steady shock 
front is recorded (a segment with AP(r) ^ const is observed directly behind the front). 
The meter type II radio emission [Gopalswamy et ai, 2009] is recorded when the forefront 
of CME 2, corresponding either to the current sheet at the frontal structure boundary or 
to the forming shock wave, is at 1.5 R© < R < 2.6 R© (crosshatched segment labelled 2 
in Figure 2). The decameter type II emission ceases to be recorded when the forefront of 
the CME, corresponding to the shock, reaches R ^ 3.5 R© (vertical line in Figure 2). 

CME 3, 3 June 2007. The difference brightness profiles AP(r) for this event are 
shown in Figure 4B. These profiles are based on the CORla data, in the frontal structure 
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coordinate system in the direction a = 0° for the two consecutive instants 09:35 (crosses) 
and 09:55 (triangles). These two profiles exhibit the evolution of the disturbed zone (shown 
by crosshatching) . There is no shock. For comparison, the dashed bold line shows the 
profile AP(r) at 09:45 plotted in the sideward direction a — —60°. The steepest segment, 
at the frontal structure boundary, in this profile has a difference brightness jump, (5P/, 
which is larger than the half of the maximum brightness on the scale 8i ~ 0.02 R©, 
comparable to the CORl spatial resolution. 

The velocity of the forefront of the disturbed zone is lower than Va (hght circles in 
Figure 2) at practically all the distances. The meter and decameter type II radio emission 
(shown by crosshatching 3 and vertical line in Figure 2) is recorded when the forefront is 
at 1.7 R0 < R< 2.6 Rq, with no shock signatures observable in the difference brightness 
profiles. 

Therefore, even though the type II radio bursts were observed in all the three events in 
question, the initial stage of the burst was not related to the shock in the CME 2 case, 
while no shock was observed in the CME 3 case. Thus, analysis of the last two CMEs - 2 
and 3 - allows one to conclude that, at least, a portion of type II bursts may occur in the 
absence of CME-excited shocks. Can the observed radio bursts be caused by a different 
source? 

6. Discussion of a possible shock-unrelated mechanism for the type II radio 

emission in the soleir corona 

Note, first of all, that a magnetic field jump ^Bi must correspond to a brightness 

jump bPi on the same spatial scale 5/ at the CME frontal structure boundary under the 
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conditions of rarefied magnetized coronal plasma. Drift current instability at the CME 
frontal structure boundary (similarly to drift current instability at the laminar shock wave 
front [Zaitsev, 1965]) may serve as the possible cause of plasma oscillations (subsequently 
transforming into radio emission). The condition Vd > Vtb must be satisfied for such 
oscillations to be excited, where and Vre are, respectively, the drift and thermal electron 
velocities [Chen, 1984]. Let us estimate their values. 

The drift velocity can be estimated from the Maxwell equation: Vd ~ SBic/AireNSi, 
where 5Bj and 6j are, respectively, the amplitude and scale of the magnetic field jump, 
c is the fight speed, N is the electron density, e is the electron charge. If the CME is 
regarded as a "magnetic barrier" moving in the solar wind, then 6i ~ u/coei in rarified 
plasma, where ouei — eSBi/c{mime)^'^ is the hybrid cyclotron frequency [Longmire, 1963]. 
Here, c is the light speed, e is the electron charge, rrii and nie are the ion and electron 
mass, respectively. Assuming u ~ Va, we obtain 5i ~ c/oUpe, where oUpe — [AirNe^ /me)^^^ 
is plasma electron frequency. This allows the electron drift velocity in the current sheet 
at the CME frontal structure boundary to be estimated as: 

Vd ^ 6Bic/A7reN6i « 6Bi/{47rNmef-^ 

Let Bq be undisturbed magnetic field directly ahead of a CME. According to the cal- 
culations by Chen [1996] and heliospheric observations by Lin et al. [2008], the ratio 
{SBj + Bq)/ Bq can be around 2-3 and, correspondingly, SBj/Bq f» 1-2. Presumably, the 
ratio {5Bi + Bq)/Bq should not vary significantly as a CME moves in the corona and 
heliosphere away from the Sun. Therefore, we obtain 8Bi ^ (1-2) Bq « (1-2) G for Bq 
1 G in the corona. 
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Taking density N ^ 3 x 10^ cm~^ as an estimate for the streamer belt at = 2 R©, 
and given SBj (1-2) G, we may estimate the drift velocity ^ (5-10) x 10^ km/s at 
the CME boundary. At the same time, the thermal electron velocity V^e ~ 5 x 10^ km/s, 
for the coronal temperature Tg !^ 1.5 x 10^ K, i.e. the condition Vd > Vxe is satisfied. 

Thus, the type II radio bursts observed in the CME 2 and CME 3 events may be caused 
by a drift current instability evolving at the CME frontal structure boundary. 
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Figure 1. A and C, slow CME 5 May 1997; B and D, fast CME 20 September 1997. A and 
B, difference brigfitness isolines, L is tlie position angle, the coordinate axes are in units of Rq. 
C and D, difference brightness distributions depending on distance r, measured from the CME 
frontal structure center (point O), along two different sections 1 and 2, whose direction is shown 
by dashed lines in the top panels. LASCO C2 data. 

Figure 2. Velocities u — V — Vsw relative to the surrounding SW depending on the distance 
from the Sun center for the CME frontal structure (light symbols) or for a shock ahead of the 
CME (black symbols) in the propagation direction. Black triangles inside the circles correspond 
to the shock front at, presumably, the forming stage. Dashed curve is velocity Vsw of quasi- 
stationary slow solar wind in the streamer belt in Wang et al. [2000]. Solid curve is the Alfven 
velocity in the streamer belt in Mann et al. [1999]. 

Figure 3. Difference polarization brightness profiles AP{R) at consecutive instants based on 
Mark 4 data for the 25 March 2008 CME. Light circle profiles correspond to the undisturbed 
corona. To the right of each instant are shown the parameters indicating the profile position in 
the image: its reference point {PAc, Rc) and direction (angle a). 

Figure 4. Difference brightness profiles AP(r) in a coordinate system tied to the frontal 
structure at consecutive instants based on CORla data: A) 31 December 2007 CME; B) 3 June 
2007 CME. The inscriptions are as in Figure 3. 
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